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1.  INTRODUCTION 


Since  publication  of  DASA-1800-II^  in  1966,  the  Defense 
Land  Fallout  Interpretative  Code  (DELFIC)  has  undergone  substantial 

/revision  in  all  of  its  modules.  These  revisions  have  created  some 

/ 

/  new  demands  on  the  Initial  Conditions  Module  (ICM)  and  removed  some 
/  i 

(  old  restrictions.  Of  most  direct  consequence  to  the  ICM  are  changes 
in  the  Cloud  Rise  Module  (CRM) ^ . 

The  new  CRM  accounts  for  wind  shear  effects  on  the  cloud 
rise  dynamics.  Therefore,  shot-time  winds  above  ground  zero  are 
input  via  the  ICM  rather  than  via  the  Cloud  Rise-Transport  Interface 
Module  (CRTIM)  as  was  done  originally.  The  old  CRM  could  accept  no 
more  than  forty  particle  size  classes,  and  the  size  class  structure 
was  rigidly  prescribed.  These  restrictions  have  been  relaxed  in  the 
new  CRM,  and  the  ICM  has  been  revised  accordingly.  In  addition,  the 
ICM  has  been  given  a  capability  to  accept  parameters  that  define  a 
power-law  particle  size  distribution  function.  From  these  parameters, 
it  constructs  a  particle  size  class  table  with  a  user-specified 
number  of  entries. 

Subroutines  LINK1  and  DSTBN  have  been  revised,  and  a  new 
subroutine,  SllWIND,  which  is  called  by  LINK1,  has  been  created.  Sub¬ 
routines  MASS,  TEMP,  TIME,  and  VAPOR  remain  unchanged. 

Subroutine  LINK1  is  tne  ICM  executive  program.  Subroutine 
DSTBN  constructs  particle  size  class  tables  for  lognormal  and  power- 
law  particle  distributions.  Subroutine  SHWIND  reads  in  the  shot-time 
winds  above  ground  zero. 

The  logic  of  the  ICM  consists  of  a  card  input,  which  is 
described  in  Table  1,  followed  by  serial  exercise  of  the  subordinate 
subroutines.  Adequate  detailed  documentation  is  provided  by  the 
FORTRAN  statement  listings. 

Use  of  the  ICM  is  quite  simple  with  one  exception:  defini 
tion  of  particle  size  distributions.  Therefore,  the  bulk  of  this 
supplement  is  devoted  to  discussions  of  particle  size  distributions. 


2.  THE  LOGNORMAL  DISTRIBUTION 


2 . 1  Fundamentals 

A  variable  x  is  said  to  be  normally  distributed  if  the 
probability  of  its  occurrt'nce  in  the  range  x  to  x  +  dx  is  given  by 


dN(x|y,o  )  =  -  exp 

a/2vT 


dx  , 


(1) 


where  y  is  the  mean  value  of  x  and  a  is  the  variance  of  x.  The 
square  root  of  the  variance,  a,  is  called  the  standard  deviation. 

To  define  a  lognormal  distribution,  we  make  the  trans¬ 
formation 


x  =  £n(y)  . 


(2) 


In  terns  of  the  variable  y  Eq.  (1)  becomes 


dA(yjy,o2)  =  1 


a/Zw 


exp 


i  &ny  -  y 1 
*2  oi 


d(Zny) 


(3) 


(3) 

and  y  is  said  to  be  lognormally  distributed. 

Some  statistical  properties  of  y  are  as  follows: 


mean(y)  =  exp(y  +  ^  t>‘) 


(4) 


median (y)  =  exp(y) 


(5) 


mode  (y)  =  exp(y  -  o  )  (6) 

variance  (y)  =  [exp (a2)  -  1  ]  exp(2y  +  a2) .  (7) 
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Le*-  and  1  he  the  geometric  mean  and  geometric  standard  deviation 
oi  y.  Then 


2_  =  median  (y)  =  exp(u)  (8) 

and 

s  =  exp(a)  .  (9) 

2 

Let  Aj  be  the  j-th  moment  of  A(y|p,o  )  about  the  origin. 
Then  by  definition 


dA(y|w,o^) 


(10) 


and  from  the  properties  of  the  normal  distribution  it  follows  that 


Aj  =  exp(ju  +  j2o2)  .  (11) 

A  feature  that  distinguishes  the  lognormal  distribution 
from  the  normal  distribution  is  the  existance  of  moment  distributions. 
The  j-th  moment  distribution  is  d* fined  as 


A(y|u,o2)j  =  yy  J  tj dA( 1 1 u,o2) 


which  can  be  shown  to  be 


(3) 


A(y| V»°2)j  =  A(y| y  +  jo2, o2)  . 


(12) 


(13) 
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A 


The  moment  distributions  provide  simple  relationships  between  log¬ 
normal  distributions  of  number,  surface  area,  and  volume  of  particles 
with  respect  to  their  diameters. 


2.2  Application  to  Particle  Distributions 

In  discussions  of  lognormal  particle  distributions,  con¬ 
fusion  frequently  arises  because  distinction  is  not  clearly  made 
between  p  and  and  between  o  and  s .  Since  particular  values  of  p 
and  o  depend  on  the  base  of  the  logarithms  used,  we  have  chosen  to 
confine  our  discussions  in  the  DELFIC  documentation  to  parameters  in 
the  form  of  and  s. 

Suppose  that  we  have  plotted  cumulative  numbers  of  parti¬ 
cles  versus  diameter  on  log-probability  graph  paper  and  have  obtained 
the  curve  shown  in  Figure  1.  This  straight-line  curve  indicates  that 
the  distribution  of  particle  number  with  respect  to  diameter,  D,  is 
lognormal.  Thus,  D  is  equivalent  to  y  in  Eq.  (3),  and  from  Eqs.  (8) 
and  (9)  we  have 


and 


S  D84.13/D50 


These  are  the  quantities  DMEAN  and  SD,  respectively,  that  are  required 
as  input  to  subroutine  LINKl,  and  that  are  printed  by  LINK1.  DMEAN 
is  expressed  in  units  of  micrometers.  (In  words,  DMEAN  is  the  median 
particle  diameter  in  the  distribution  of  numbers  of  particles  with 
respect  to  their  diameters.)  If  the  user  specifies  a  lognormal  dis¬ 
tribution  but  does  not  input  values  for  -DMEAN  and  SD,  the  program 
supplies  the  values 
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DMEAN  =  x  =  0.407  pm 
SD  =  s  =  4.0 


As  noted  above,  the  properties  of  the  moment  distribu¬ 
tions  are  useful  in  interrelating  distributions  of  particle  number, 
surface  area,  and  volume  with  respect  to  particle  diameter.  This  is 
because  the  number  distribution  is  the  zeroth  moment  distribution 
with  respect  to  diameter,  surface  area  is  distributed  via  the  second 
moment  distribution,  and  volume  is  distributed  via  the  third  moment 
distribution.  Thus,  if  we  assume  spherical  particles  and  if  the 
parameters  p  and  a  are  known  for  either  the  particle  number,  or  parti¬ 
cle  area,  or  particle  volume  distribution  with  diameter,  then  the 
other  distributions  can  be  determined  from  the  equations  below.  The 
parameter  o  is  the  same  for  all  three  distribution!?.  If  we  use  N,  S, 
and  V  as  subscripts  to  denote  number,  surface  area,  and  volume, 
respectively,  we  have  fron  Eq.  (13) 


ps  =  uN  +  20 


PV  =  WN  +  30 


where  p  and  o  are  related  to  %  and  s  by  Eqs.  (8)  and  (9). 

If  base  10  logarithms  are  used  instead  of  natural  loga¬ 
rithms,  we  distinguish  the  distribution  parameters  by  use  of  primes, 
p'  and  o',  and  the  relations  become 

p’  -  PjJ  +  2M10)(o')2 

+  p^  +  3M10)(o')2 


where  «.n(10)  =  2.3026. 
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The  distribution  of  particle  mass  with  respect  to  diameter 
is  taken  to  be  equivalent  to  the  volume  distribution  with  respect  to 
diameter.  This  implies  that  all  particles  have  the  same  density. 

2.3  Particle  Size  Class  Tables 

For  computation  purposes,  the  continuous  lognormal  distri¬ 
bution  is  replaced  by  a  histogram.  The  computer  program,  via  sub¬ 
routine  DSTBN,  does  this  automatically  by  use  of  the  distribution 
parameters  and  the  number  of  size  classes,  NDSTR,  which  is  input  by 
the  user. 

The  user  specifies  parameters  DMEAN,  SD,  and  NDSTR.  From 
these,  the  parameters  u  ,  o,  and  are  determined  via 

UN  =  £n( DMEAN) 

o  =  £n(SD) 

=  yN  +  3°2  * 

Define  the  normal  distribution  function  argument  x  as 

£n(D)  -  uv 

x  =  - 

o 

where  D  is  particle  diametei .  Then 

x 

N(x)  =  — —  f  exp(-t  /2)dt 

/2tT  J 
—00 

Subroutine  DSTBN  constructs  the  p&.  cle  size  class  table 
(i.e.,  histogram)  as  follows.  Each  size  class  contains  a  constant 
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volume  fraction,  AN^,  of 

ANy  =  1/NDSTR  . 

Let  D^,  i  =  1,2 . ^NDSTR,  be  the  upper  (i.e.  the  larger  particle) 

boundary  diameter  of  the  i-th  particle  size  class.  The  table  is 
ordered  with  the  largest  particles  in  the  first  size  class,  and  so  on. 
Then,  for  the  i-th  size  class 

N(x.)  =  iAN 
1  V 

and 

^i+P  =  V  +  uv  ' 

The  upper  boundary  of  the  first  size  class,  D^,  and  the 
lower  boundary  of  the  last  size  class,  DNnSTR+1>  are  special  cases. 
These  are  taken  to  be  the  diameters  at  AN^/2  and  l-ANy/2,  respectively. 
That  is, 


and 


N(xx) 


n(xndstr+i->  1 


In  these  calculations  x  is  determined  from  given  N(x) 
via  equation  26.2.23  of  Reference  4. 

The  central  particle  diameter  for  the  i-th  class,  d^,  is 

given  by 


d. 

x 


"  V  DiDi+l  ' 


I 

’ft 
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i 


If  NDSTR.  =  1,  a  single  size  class  is  created  with 

•  *  i 

Dx  =  (DMEAN)  *  (5  J  *  SD)  , 

^  \ 

D2  =  (DMEAN)  /  (5.0  *  SD) 

;  I 

i  , 

and  ;  i 


d1  =  DMEAN  . 


i 


l  ‘ 
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I 


! 
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» — ^  - 
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3.  THE  POWER-LAW  DISTRIBUTION 


3.1  Fundamentals 

Mathematically  speaking,  power-law  distributions  are 
meaningless  since  distribution  functions  cannot  be  defined  for  them. 
This  is  because  the  power-law  function  is  not  properly  bounded  for 
zero  argument  Freiling  has  shown  that  fallout  particle  distribu¬ 
tions  that  have  been  represented  by  power-law  functions  can  equally 
well  be  fitted  by  lognormal  distribution  functions. ^  The  implica¬ 
tion  of  Freiling's  work  is  that  power-law  distributions  would  be  more 
accurately  described  as  truncated  lognormal  distributions.  Neverthe¬ 
less,  power-law  distributions  frequently  are  useful  in  fallout  work. 
Define  the  power-law  frequency  as 

] 

df(D|k,X)  =  kD~XdD  ,  (14) 

where  df(D|k,X)  is  the  number  of  particles  in  the  diameter  range  D  to 
D  +  dD.  If  we  assume  spherical  particles  with  constant  density,  p, 
we  have 


dF 


irpk 
6M  ’X 


i  n3_x 

TTpkP 

6M 


dD  , 


(15) 


where  dF|D  |^,x|  is  the  fraction  of  the  total  fallout  mass,  M,  in 
the  diameter  range  D  to  D  +  dD. 

The  mass  fraction  of  particles  in  the  macro-range  from  Dt 
to  0^  is  obtained  by  integration  of  Eq.  (15)  between  these  limits  to 
give 


A-p  -  ”Pk 

'  6M(4-X) 


0  <  X  <  4 


(16) 
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3.2  Particle  Data  Analysis 

Suppose  that  we  have  obtained  a  sample  of  fallout  particles. 
We  weigh  the  sample  to  obtain  M  (kg),  and  we  size  the  sample  into  N 
fractions,  the  i-th  fraction  containing  particles  in  the  diameter 
range  AD.  centered  on  D^  (meters).  We  weigh  each  fraction  and  obtain 

the  mass  fractions  AF..  We  determine  that  the  average  particle 
3  1 

density  is  p  (kg/m  ) . 

To  obtain  the  power  law  distribution  parameters  k  and  X, 
we  plot  log(AFi/AD^)  versus  log(D^) .  A  straight  line  is  fitted  to  the 
data.  From  Eq.  (15),  we  see  that  the  intercept  and  slope  are 

intercept  =  c  =  logj^pj  , 

and 


slope  =  m  =  3  -  X 


Then 


and 


X  =  3  -  m 


,  6M  .  -1.  . 

k  =  —  Tog  (c) 


When  X  and  k  are  determined  from  M  expressed  in  kilograms, 

3 

D  and  AD  in  meters,  and  p  in  kg/m  ,  they  can  be  input  to  subroutine 
LINK1  as  EXPO  and  CAY,  respectively. 

3.3  Particle  Size  Class  Tables 

For  use  in  fallout  calculations,  subroutine  DSTBN  creates 
a  histogram  representation  of  the  power  law  distribution.  The  histo¬ 
gram  is  comprised  of  NDSTR  particle  size  classes,  where  NDSTR  is 
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specified  by  the  user.  The  mass  fraction  in  each  size  class,  AF,  is 
the  constant 

AF  =  1/NDSTR  . 


Let  be  the  upper  (i.e.  larger  particle)  boundary  of  the 
i-th  particle  size  class.  The  table  is  ordered  with  the  largest 
particles  in  the  first  class,  and  so  on.  Then  the  smallest  particles 
are  contained  in  the  NDSTRth  class.  If  we  assume  that  the  smallest 
particle  in  this  class  is  much  smaller  than  we  see  from 

Eq.  (16)  that 


4-X  _  6M(4-X) 

UNDSTR  upk 


AF  . 


By  recursive  use  of  this  relation  with  Eq.  (16),  we  find  that 

1 

Dt  -  (NDSTR  -  i+l)4-X  Dbdstr  . 


Size  class  central  diameters,  d^,  are 


di  -V“A 


i+1 


To  establish  a  central  and  lower  boundary  diameter  for  the 
NDSTRth  class,  we  say  that 

.  _  M2^  n 

dNDSTR  \ 2 }  °NDSTR 
and 


°NDSTR+1  ~  (dNDSTR^  /DNDSTR  * 
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4.  TABULAR  DISTRIBUTIONS 


4.1  Particle  Size  Class  Tables 

If  the  user  so  desires,  he  can  input  his  particle  size 
distribution  in  histogram  form  with  NDSTR  size  classes.  The  table 
of  size  classes  must  be  arranged  in  descending  order  of  particle 
diameter.  Each  size  class  is  defined  in  the  input  by  its  upper 
(i.e.  larger  particle)  boundary  diameter,  D^,  and  mass  fraction,  AF^ 
These  two  data  are  punched  on  a  separate  card  for  each  size  class. 

The  last  card  in  the  deck  contains  the  lower  boundary  diameter  of  the 
NDSTRth  size  class.  Central  particle  diameters,  d^,  are  computed  to  be 

di=  (Di  +  IW/2  • 


Preceding  page  blank 


5.  USER  INFORMATION 


5 . 1  Card  Input 

The  I CM  card  input  is  described  in  Table  I.  This  table 
and  the  discussions  in  Sections  2.2,  2.3,  3.2,  and  4.1  provide 
adequate  information  for  use  of  the  code. 

5.2  Output 

Though  the  printed  output  has  been  modified  somewhat, 
the  example  output  presented  in  DASA-1800-II  is  still  satisfactory. 

Communication  with  the  Cloud  Rise  Module  is  via 
C0MM0N/SET1/.  The  contents  of  C0MM0N/SET1/  is  described  in  Table  2.2 
of  DASA-1800-III  (Revised). 


Preceding  bm 
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TABLE  I 

CARD  INPUTS  TO  THE  INITIAL  CONDITIONS  MODULE 


Card  Number 

Concents 

Variable  Names 
and  FORMATS 

1 

ICM  Run  Identifier 

DETID(J) ,  J=l,12 
(12A6) 

2 

Control  integer  to  specify  particle  size 
distribution  type: 

1  lognormal 

2  power -law 

3  tabular 

IDISTR  (15) 

3 

Number  of  particle  size  classes 

NDSTR  (15) 

4(a)* 

(For  lognormal  particle  size  distribution) 
Explosion  yield  (KT) ,  height  of  burst 
above  GZ(m),  soil  class  indicator: 

1.0  for  siliceous 

2.0  for  calcareous, 

median  particle  diameter  (pm) ,  geometric 
standard  deviation  of  the  particle  size 
distribution,  and  particle  density 
(g/cm3).  (See  Section  2.2.) 

W,  HEIGHT,  USOIL, 
DMEAN,  SD,  DNS 
(6F10.3) 

4(b)* 

(For  power-law  particle  size  distribution) 
Yield  (KT) ,  height  of  burst  (m) ,  soil 
class  indicator  (see  above),  exponent  in 
the  particle  size  distribution  frequency 
function,  coefficient  in  the  particle  size 
distribution  frequency  function,  particle 
density  (g/cm^).  (See  Section  3.2.) 

W,  HEIGHT,  USOIL, 
EXPO,  CAY,  DNS 
(6F10.3) 

4(c)* 

(For  a  tabular  particle  size  distribution) 
Yield  (KT) ,  height  of  burst  (m) ,  soil  class 
indicator  (see  above),  particle  density 
(g/cm3) . 

W,  HEIGHT,  USOIL, 
DNS  (4F10.3) 

,  ** 
4(c)1 

• 

• 

A  table  of  upper  boundary  particle 

diameters  (pm)  and  mass  fractions 

• 

DIAM(J) ,  FMASS(J) , 
J=l, NDSTR  (2E12.5) 
• 

• 

** 

4(c)ndstr+i 

• 

The  lower  boundary  diameter  (pm)  of  the 
last  particle  size  class.  (See  Section  4.1.) 

* 

DIAM(NDSTR+1) 

(E12.5) 
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Table  I  (continued) 


Variable  Names 

Card  Number _ Contents _ and  FORMATS 

5  Number  of  entires  in  the  wind  data  table  NHODO  (15) 

icick 

6  For  each  entry  in  the  wind  data  table:  ZV(J),  VX(J), 

altitude  (m,  relative  to  msl),  x  com-  VY(J),  J=l>  NHODO 

ponent  of  wind  (m/sec),  y  component  of  (3F12.3) 

wind  (m/sec) 


*  One  of  the  cards  4(a),  4(b),  or  4(c)  is  read  according  to  whether 
IDISTR  is  1,  2,  or  3. 

**  These  cards  are  read  only  for  a  tabular  distribution. 

***  These  cards  are  read  only  if  NHODO  >  0. 
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6.  FORTRAN  STATEMENT  LISTINGS 


Complete  FORTRAN  statement  listings  are  given  for  the 
following  subroutines.  These  subroutines  are  operational  on  the 
UNIVAC  1108. 


SUBROUTINE 

Page 

LINK1 

22 

DSTBN 

27 

SHfIND 

29 

The  machine  used  to  prepare  these  listings  prints  a 
//  symbol  to  represent  a  4-8  punch;  this  symbol  should  be  an 
apostrophe  (').  In  FORMAT  and  DATA  statements,  the  apostrophe  is 
used  to  define  Hollerith  character  fields. 


Preceding  page  blank 


21 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


SUBROUTINE  LTNKI  LTNKI 

INITIAL  CONDITIONS  (FIREBALL)  H00UL-.  LINKI 

NT.  AUBURN  RESEARCH  ASSOCIATE^  jANUARf  1972  LINKi 

LINKI 

I.TNK1 

PROGRAM  TO  DETERMINE  THE  INITIAL  CONDITIONS  SPECIFICATIONS  OF  LTNKI 

TIME,  TEMPERATURF,  TOTAL  SOIL  MASS*  FRACTION  OF  THF  SOIL  aURDFN  INLINKI 
THE  VAPOR  PHASE,  ANC  THE  SI7E  FREQUENCY  DISTRIBUTION  OF  THF  LINKI 

CONDENSfcO  PHASE  SOU  LTNKI 

LINKI 

THF  FIRST  CARD  CONTAINS  ANY  ARBITRARY  ALPHANUMERIC  IDENTIFICATION. LINKI 

LINKI 


OTHER  INPUT  PARAMETERS  ARE  -  TECT  PARAMETER  UDTSTR)  TO  DETERMINE  LINK) 
IF  THE  PARTICLE  SIZE  RRFOUENCY  niSTRIBUTICN  IS  LOG-NORMAL,  POWER  LINKI 
LAW,  OR  TABULAR,  Y I  FLO  IN  KILOTONS,  HE IGHT (DEPTH )  OF  BURST  IN  LINKI 
METERS,  A  SOIL  TYPE  INDICATOR,  FAuLOUT  PARTICLE  DENSITY ( GM/CM**3 ) , LINK! 
MEANIHICPOMETERS)  AND  STANDARD  DEVIATION  FOP  A  LOG-NO»MAL  PART IOLELTNK1 
SIZE  FREQUENCY  DISTRIBUTION,  THF  NUMBER  OF  PARTICLE  SIZE  CLASSES  LINKI 
IN  THE  PARTICLE  SIZE  FREQUENCY  DISTRIBUTION.  IF  EITHER  A  TABULAR  LINKI 


OR  POWER  LAW  DISTRIBUTION  IS  USED  ,  THE  MEAN  AND  STANDARD  LINKI 
DEVIATION  ARE  NOT  CCf.LFP  FOP  SINCE  THEY  DC  NOT  APPLY.  IF  A  LINKI 
LOG-NORMAL  DISTRIBUTION  IS  TO  BE  SUPPLIED  BY  THE  PROGRAM,  THE  LINKI 
MEAN  AND  STANDARD  DEVIATION  FIELDS  ARE  LEFT  BLANK.  LINKI 
SHOT  TIME  WINDS  ABOVE  G7  ALSO  ARE  INPUT.  THESE  ARE  USED  TO  LINKI 
COMPUTE  WIND  SHEAR  EFFECTS  ON  CLOUD  RISE  AND  FALLOUT  ADVECTION  LINKI 
DURING  THE  CLOUD  RISE  TIME  INTERVAL.  LINKi 


LTNKi 

FOR  UNDERGROUND  BURSTS  INPUT  DEPTH  OF  BURST  AS  A  NEGATIVE  NUHREP  LINKI 

LINKI 

THE  OUTPUT  UNITS  ARE  MASS  IN  KILOGRAMS,  LENGTH  IN  METERS,  TTME  IN  LINKI 


SECONDS,  TEMPERATURE  IN  DEGREES  KELVIN,  YIELO  IN  KILOTONS,  LTNKI 

DISTRIBUTION  PARAMETERS  IN  MICRONS  LTNKI 

LINKI 

*•**«»»»»»»«» »**»*•»  GLOSSARY  *********************** *********** LINK! 

LINKi 

CAY  COEFFICIENT  OF  THE  FREQUENCY  FUNCTION  FOR  THE  POWFR  LINKI 

LAW  PARTICLE  SIZE  FREQUENCY  DISTRIBUTION  LINKI 

DETIDm  INITIAL  CONDITIONS  IDENTIFICATION  ARRAY  LINKI 


OIAM(I)  ARRAY (201) ,  UPPER  BOUNDARY  OF  THE  I-TH  PARTICLE  SIZE  LINK  1 

CLASS.  THE  LAST  ENTRY  IN  THE  DIAM  ARRAY  IS  THE  LOWERLINK1 


BOUNDARY  CF  THE  LAST ! SMALLEST)  PARTICLE  SIZE  CLASS.  LTNKI 

THE  LENGTH  OF  THE  OIAM  ARRAY  IS  ALWAYS  ONE  GREATER  LIN<1 

THAN  THE  NUMBER  OF  SIZE  CLASSESlMICPOMETf RS)  LTNKI 

OMEAN  MEDIAN  DIAMETER  (MICROMETERS)  OF  LOGNORMAL  P ARTICLF  LINKI 

SIZE  DISTRIBUTION  LINKi 

QMS  FALLOUT  PARTICLE  DENSITY  (GM/CM**3)  LTNKI 

EXPO  EXPONENT  CF  THE  FREQUENCY  FUNCTION  FOR  THE  POWFR  LINK! 

LAW  PARTICLE  SIZE  FREQUENCY  DISTRIBUTION  LINKI 

FMASS(I)  ARRAY  OF  FRACTION  OF  TOTAL  PARTICULATE  MASS  IN  I-TH  LINKI 

PARTICLE  SIZE  CLASS.  MAXIMUM  LENGTH  OF  ARRAY  =  200  LINKI 
HEIGHT  HEIGHT  OF  BURST  <HETERS>  ABOVE  GROUND  ZERO  LINKI 

IOISTR  CONTROL  INTEGER  FOR  PARTICLE  SIZE  DISTRIBUTION  LTNKi 

1  -  LOGNORMAL  DISTRIBUTION  LINKI 

2  -  POWER  LAW  DISTRIBUTION  LINKI 

3  -  TABULAR  DISTRIBUTION  REAO  IN  ON  CAROS  C ARRAY  WHY) LINKI 

IS  CONTPOL  INTEGER  SPECIFIES  WHETHER  LOGNORMAL  LINKi 

DISTRIBUTION  IS  SPECUIED  BY  THE  USER  OR  BY  THF  LINKI 


1 

? 

3 

4 

0 

6 

7 

5 
a 

10 

11 

1? 

I? 

14 

lc 

16 

17 

IP 

IP 

20 

21 

22 

23 

2A 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
3E 

37 

38 

39 

40 

41 
4? 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 
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o  o  o  o  oooooooooonoooonoooooooociooooooo 


I S  In 

rsnuT 

NDSTR 

NH000 

PS(I) 


the 

TMPl 

TMP? 

r?M 

USOTL 


VPR 
V*(  I) 

V  Y  ( I ) 

W 

ZSCl 

ZV(I) 


PROGRAM 

0  -  PROGRAM  S°EC IFIEO  LOG-NORMAL  DISTRIBUTION 
1  -  USER  SPECIFIED  LOG-NORMAL  DISTRIBUTION 
SYSTEM  INPUT  TAPE 
SYSTEM  OUTPUT  TAPE 

NUMBER  OF  PARTICLE  SIZE  CLASSES.  MAXIMUM  =  200 
NUMRE 0  CP  ENTRIES  TN  WIND  HODOGRAPH  TABLE 
ARRAY  (  200).  PARTICLE  SIZE  CLASS  MIDPOINT  01 AMETFR 
(METERS) 

STANOARO  DEVIATION  OF  LOGNORMAL  PARTICLE  SIZE 
D I STRI CUT I ON (DIMENSIONLESS) 

MASS  OF  CONDENSED  PHASE  MATERIAL  AT  SPECIFICATION 
TIME 

TIME  OF  INITIAL  CONDITIONS  SPECIFIC AT  ION 
AVERAGE  TEMPERATURE  OF  GAS  IN  ClOUD 
AVERAGE  TEMPERATURE  OF  CONDENSED  PHASE  MATF°IAL  IN 
CLOUD 

TEMPORARY  STORAGE 
SOIL  CLASS  INDICATOR 
1.0  FOP  SILICEOUS 
2.0  EOF  CALCAREOUS 

MASS  OF  VAPOR  IN  CLOUO  AT  SPECIFICATION  TIME 
X-COMP CNENT  OF  WINO  VELOCITY  AT  MIND  HODOGPAPH 
STRATUM  I  (METERS/SEC.) 

Y-COMRCNENT  OF  HIND  VELOCITY  AT  WIND  HODOGRAPH 
STRATUM  I  (METERS/SEC.)  * 

WEAPON  YlELn  (<T) 

SCALED  HEIGHT  OF  BUPST 

ALTITUCE  CF  THE  WIND  COMPONENTS  VX<1)  AND  VY(I) 


VX( 1)  AND  V Y (I) 


(METERS  RELATIVE 


MFAN  SEA  LEVEL) 


COMMON  /SET  1/ 

1CAY  ,DETID(12) 

2FMASS (200) , IOISTR  , IEXEC  ,IRISE  ,IStN  ,ISOUT  ,LTNK1  93 

3NOSTP  r  PS ( 2  0  0 )  ,SO  ,SSAM  ,TME  ,TMP1  , LTNK1  94 

4TMP2  ,T2M  ,USOIL  ,VPR  ,W  , HEIGHT  ,LINK1  95 

5ZSCL  ,NHODO  t IV (200 )  ,VX(200)  ,VY(200)  LIN<1  96 

LINK1  97 
LIN<1  98 

LTN<1100 

F0PMAT(12A6)  LINK1101 

’  FORMAT (/3X»60HTHE  SPECIFIED  STANDARO  DEVIATION  IS  NEGATIVE  HENCE  ILINK1102 

INCORRECT///)  LINK1103 

1  FOPMAT (7F10 .3)  LINKU04 

>  FORMAT (///25X28H****  INPUT  PARAMETERS  * **V?0X , 5HYIELD ,4QX ,E12.5LIN<1105 

1»2X,2HKT/20X»24H HEIGHT  OP  DEPTH  0*  BUR;'T,21X,E12.5,2X,6HMETE  RS/20XLINK1106 
2 » 13HS0IL  CATEGORY)  LTNK11Q7 

FORMAT (1H*  ,65X  , 9HS ILICEOUS )  LINK1108 

»  FORHATdHM,  65X,10HCALCAPEOUS)  LTNK1109 

FORMAT (//20X,  36HPART ICLE  SIZE  FREQUENCY  DISTRIBUTlON/LINidllO 

125X32HA  LOG-NORMAL  CISTRIBUTION  WITH  - /3QX , 15HME0IAN  DIAMETER, 20  ,LINK1111 
2E12.5, 2X, 11 HMICROMFTERS/?OX,29H GEOMETRIC  STANDARD  DEVIATION,  7X,  LTN<1112 

3E12.5/25X,  34HTHIS  DISTRIBUTION  HAS  SPECIFIED  BY)  LTNK1113 

i  FORMATI1H* ,65X,11HTEE  PPOGRAM)  LINK1114 

I  F0RMAT(1H+,65X,8HTHE  USER)  LlN«dU5 

.0  FORMAT  (15)  LTNKU16 


,CTAM(201) 
, IEXEC 
,SO 

,USOIL 

t ZV (200 ) 


, DMt AN 
, IRISE 
,  S  S  AM 
,  VPR 
» VX  (200) 


,DNS 

,ISIN 

,TME 

»W 

, VY  (2  00) 


,EXPO 
, ISOUT 
,TMP1 
♦HEIGHT 


LINK1  59 
LTNK1  60 
LIN<1  61 
LTNK1  62 
LINrfl  63 
LTNK1  64 
LIN<1  65 
LINK1  66 
LTN<1  67 
LINK1  68 
LIN*1  69 
LIN<1  70 
LINK1  7 1 
LTN<1  72 
LTNK1  73 
LTNK1  74 
LIN<1  7C 
LTNK1  76 
LTN<1  77 
LINK1  78 
LINK1  79 
LIN<1  80 
LINK1  81 
LINi<1  82 
LINK1  83 
LTNKI  84 
LINK1  85 
LINK1  96 
LINK1  87 
LTNKI  88 
LlNKl  89 
♦LINK1  90 
LTN<1  91 
, LINK1  92 


23 


oo  no  oonoo 


11  FORMAT  (/3X,53HTHE  SCALED  DEPTH  OF  BURST  IS  BEYOND  THE  SCOPE  OF  THFLTNK1117 

1  MODEL///)  LTNK111  3 

12  FORMAT (/3X,111HTHF  SCALED  HEIGHT  OF  BURST  IS  SUCH  THAT  THERE  IS  NOLlNKlll® 

1  SOIL  PASS  FNTRAINFC  IN  THE  CL0U0  AND  HENCE  NO  LOCAL  FALLOUT///)  LTNK1120 

13  F0PMAT1///25X37H****  INITIAL  CLOUD  PROPERTIES  AT  H  «-E12.5,14H  SFCLTN<1121 

1 CMOS  ****//20X, 23HAVERAGF  GAS  TEMPERA TURE3RX , El  2.5, 2X, 1 4HDFGRFE S  LINK112? 
2KELVTN//20X,56HAVERAGE  TEMPERATURE  OF  CONDENSED  PHASE  MATERIAL  IN  LTNK112? 
3CLOUO,5X,E12. 5, 2X,14H DEGREES  <ELVIN.'/20X,31HMASS  OF  VAPORIZED  SO ILLINK1124 
4  IN  CLCUD,30X,E12. 5 , 2X , 9HK TLOGP AHS //2 0X41HMASS  OF  CONDFNSEO  PHASE  LTNK1125 
5MATFRIAL  IN  CLOUO , ? OX , F 12. 5, 2X ,9H< ILOGR AMS)  LINK1126 

15  FOPMAK1X,  14HLEAVINT  LINK  1)  LINK1127 

16  FOPMAT (A1A///51X, A*  ********  *///l2X,AT  HF  DEPART  NLTNK112B 

1FNT  OF  DEFENSE  FALLOUT  PRFDICTIO  NLINK1129 

2  S  Y  S  T  F  MA//51X,/*  ********  *<////43X,ElNITIAL  ooNOITILINKHTO 

30NO  (FIREBALL)  MODULE A,/// 55 X , /ORE P APED  BY//46X,/MT.  AUBURN  PESE APLINK1131 
4CH  ASSOCIATES  A  755 X ,  /NEWTON,  MASS. A////25X,/****  INITIAL  CONDITIONLIN<1132 
5S  IDENTIFICATION  * ***//25X, 12A6)  LINK1133 

17  FORMAT (/3X,60HTHF  SFECIFIEO  MEAN  PARTICLE  SI7E  IS  NEGATIVE  HENCE  TLTNK1134 
INCORRECT///)  LINK113* 


18  FORMAT <//20X, 36HPARTICLF  SIZE  FREOUFNCY  DISTRIBUTION/  LINK113F 

125X41H A  TABULATED  EMPIRICAL  DISTRIBUTION  WITH  -/ 30X12 , 2X , 21 HPAPT ICLINK1 13 7 
2LF  SIZE  CLASSES)  LINK1136 

192  F0RMAT(//51X,19H*  ********  *// »  LINK1139 

193  FORMAT(/l/,9X, APAPTICLF  SIZE,  LOWER  SIZE  INTERVAL  BOUNDARY, HASS  FRLINK1140 
1EQUFNCY,  AND  UPPER  SI7F  INTERVAL  BOUNOARY//10 X , /FOR  USE  IN  PSXP, TRLTNK1141 
2ANSP0RT,  AND  ACTIVITY  CALCULATIONS  (01 AMETERS  IN  METERS) A//20X ,/ DTLINK114? 
3AMFTER<j4X, /LOWER  BfUNOARY A, 3X , /MA SS  FRACTION/ , 3X , /UPPER  BOUNnA® YALINK1143 
4//)  LINK 11 4 4 

197  FOPMAT(//20X, /PARTICLE  SIZE  FREQUENCY  DISTRIBUTI0N//25X, /POWER  L AWLINK1145 
1  DISTRIBUTION  WITH  -  A/3 OX , I 3 , IX, / P ART ICLE  SI7£  CLASSES// 30 X , /THF  LINK1146 
2SPFCIF TFD  PARAMETERS  ARF//30 X, /CAY  =/,2X,E12.5/30X,/EXPO  =/, 2X,E12LINK1147 


3.5) 

194  FORMAT (12X,I3,4(3X,F12.5)) 
195  FOPMAT (2E12 .5) 

198  FORMAT (/3X,58HTHE  PARTICLE 


SIZE  DISTRIBUTION  TABLE  IS  IMPROPERLY 


LINK114B 

LTNKH49 

LTNK1150 

OLINK1151 


2RDERED///) 

LINK115? 

LINK1153 

»** 

****LTN<1154 

*** 

*+ **LI NK1 155 

LTNK115F 

READ  INITIAL  CONDITIONS  RUN  IDENTIFIER 

LTNK1157 

READ  (ISIN,1) (OETIOCJ) , J=l,12) 

LTNK115  3 

LINK1159 

WRITE  OVERALL  TITLE 

LINK1160 

WRITE  (ISOUT, 16) ( DE TID ( J) , J= 1 , 12) 

LINK1161 

20 

REAO  (ISIN,10) IDISTR 

LINK1162 

REAOdSIN,  10)  NOSTR 

LTNK1163 

IF(NDSTR)40 1,4 01,40? 

LINK1164 

401 

NDSTP=100 

LINK1165 

402 

GO  TO  (210,220 ,211)  , IDISTR 

LINK1166 

210 

READ(ISIN,3)W,  HEIGH T,UFOIL,DME AN, SD, DNS 

LINK1167 

WAS  A  PRESHOT  PARTICLE  LOG-NORMAL  DISTRIBUTION 

SPECIFIED  BY 

LINK1168 

THE  USER  YES  TO  22 

LINK1169 

IF (0MEAN)21,21,22 

LINK1170 

21 

IS=0 

LINK117 1 

GO  TO  23 

LINK1172 

22 

IS  =  1 

LINK1173 

GO  TC  ?3 

LINK1174 

24 


oooonoo  ooo  oo 


22  0 
I'll 


C 

C 


21*5 

C 

C  23 
23 
C 

60 


63 

66 


90 

91 


C 

C 


301 

302 
305 
309 


PEAP(  ISTN,3)W,HF.IGHT,USPTL,EXP0,CAY,0NS 

LTNK1175 

GO  TO  23 

LINK1176 

REA0(ISIN,  31 W,  HETGHT,USOIl,nNS 

LINK11 77 

RE  AO  ( ISIN, 1 95 ) (DIAM(I) ,FMASS(I) ,I=1,NDSTRI 

LINK1178 

LD=NDS  TP  + 1 

LINK1179 

REA9(ISIN,195)DIAM(tD> 

LTNK11AC 

LTNKim 

CHECK  ORDERING  OF  THE  HISTOGRAM  TABLE 

LINK116? 

DO  215  1=2, ID 

LINK1183 

IFCDIAMd)  .LT.  OI AM(T-l) )  GO  TO  215 

LTNK1184 

WRITE (  ISOUT, 198) 

LINKH85 

GO  TO  190 

LINKH86 

CONTINUE 

LV-NKlt  17 

LTNY1186 

CONVERT  HOO  -  OOD  FROM  MFTERS  TO  FEFT 

LTNK1169 

HE IGHT=HEIGHT /  0 . 30 4  F 

LINK1190 

7SCL  IS  THE  SCALED  FOP  -  OOP 

LINK1191 

ZSCL  =  HEIGHT/’UW)  **  (1.073.4) ) 

LINK1192 

LINK1193 

TEST  THE  DATA  TO  SEE  IF  THE  MODEL  IS  APPROPRIATE 

LINK1194 

IF (HEIGHT) 66,66,63 

LINK1195 

IF (7SCI-190.Q) 70,70,150 

LTNK1196 

IF (ZSCLH20.0) 143,70  ,70 

LTNK1197 

CALL  TIME 

LINK1198 

CALL  TEMP 

LINK119° 

CALL  MASS 

LINK1200 

CALL  VAPOR 

LINK1201 

GO  TO  (90,95,99) , IOISTP 

LINK1202 

LINK1203 

TEST  FOR  ACCEPTABLE  SPECIFICATIONS  OF  PRE-SHOT  PARTICLE 

SIZE 

LINK1204 

FREQUENCY  DISTRIBUTION. 

LINK1205 

IF (SO) 91,92,9? 

LTNK1206 

WRITE  (ISOUT,?) 

LINK1207 

GO  TO  190 

LTNK120  3 

IF(DMEAN)94,95,95 

LINK1209 

WRITE  (ISOUT, 17) 

LINK1210 

GO  TO  190 

LTNK1211 

LINK1212 

CALL  DSTBN 

LTNK1213 

LINK121 4 

CONVERT  HOB  -  DOB  PACK  TO  METERS  FROM  FEET 

LINK1215 

HE IGHT  =  HEIGHT*  0. 30  4  6 

LTNK1216 

LINK1217 

CONVERT  VPP  AND  SSAM  cROM  GRAMS  TO  KILOGRAMS 

LINK121 3 

VP R=VPR/1000.0 

LINK1219 

DURING  COMPUTATION  SSAM  CONTAINS  THE  VALUE  OF  THE  TOTAL 

MASS  OF 

LINK1220 

GAS  AND  CONDENSED  PHASE  MATERIAL  IN  THE  CLOUD. 

LINK1221 

ssam=ssam/iooo.o-vpf 

LINK1232 

LINK1223 

WRITE  INITIAL  CONDITIONS  RESULT*? 

LINK1224 

WR TTE  (ISOU T, 4) W, HEIGHT 

LINK1225 

IF  (USOIL-1  •  0)  301,30  .1,302 

LINK1226 

WRITE  (ISOUT, 5) 

LINK1227 

GO  TO  305 

LINK122B 

WRITE  (ISOUT, 6) 

LINK1229 

GO  TO  (309,310,311) ,IOISTR 

LINK1230 

WRITE (ISOUT, f )ONE AN,  SO 

LINK1231 

IF  (IS)  102, 10  3, 102 

LINK1232 
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r>  o  o 


103  WRITE  (ISOUT, M 
GO  TO  315 

10?  WRITE  (ISOUT, 9) 

GO  TO  315 

311  WRITE(TSOUT,13l MOST  R 

PRINT  FINAL  PARTICLF  SIZF  CLASS 

315  WRITE (IS OUT, 193) 

00  60?  J=1 , NOSTR 
JO=J+l 

0Ml=0IAM(J0)*1.0E-6 
OH2=niAM(  J)  n.  0E-6 

602  WRITE (ISOUT, 194) J,PS(J) ,0M1, FMASS ( J) , DM2 
GO  TO  106 

310  WRITE (ISOUT, 197) NOSTR, CAY, EXPO 
GO  TO  315 

106  WRITE(ISOUT,13) THE , TWP1 , TWP2, VPR, SSAM 

113  WRITE ( ISOUT ,192) 

200  CALL  SHWINO 

WRITE (ISOUT ,151 
RETURN 

143  WRITE  (ISOUT, li» 

GO  TO  190 

150  WRITE  (ISOUT, 121 

190  CALL  EXIT 
ENn 


LTNK1233 

LINK1234 

LIN*1235 

LTNK1236 

LTNK1237 

LINK1233 

LINK1239 

LIN<1?40 

LTNK1241 

LTNX124? 

LINK1243 

LINX1244 

LINX1245 

LINK1246 

LTNK1R47 

LINK1248 

LIN<124° 

LINK1250 

LTNX1251 

LTNK125? 

LINK1253 

LINK1254 

LINK1255 

LINK1256 

LINK1257 

LINi<1258 

LINX125R 


ooo  oooo  oooooooo 


I 


SUBPOUTINE  OSTBN 
COMMON  /SET 1/ 


1CAY 

,0ETID<12) 

, OIAM  <201) 

, DMEAN 

,DNS 

,FXPO 

?FMASS<200) 

,  IOISTR 

, IEXEC 

, irise 

,  IS  IN 

,  TSOUT 

3NDSTR 

, PS <200) 

,  SO 

,ssam 

,TME 

,  TMP1, 

4TMP2 

,  T2M 

,USOIL ‘ 

,VPR 

'  ,w 

, HEIGHT 

5ZSCL 

,NHOnn 

,?V<?00) 

,  VX  <  2  0  0 ) 

,V  Y  <2  00) 

! 

LOGNORMAL  DISTRIBUTION  TO  100 

POWER  FUNCTION  DIS  TFIBUTION' TO  200  . 

TABULAR  OtSTRIBUTIOM  TO  700 

EOIJATI0N  26.2.23  OF  NBS-AMS  55  HANDBOOK  IS  USED  TO  COMPUTE  THE 
PROBABILITY  FUNCTION  ARGUMENT  FROM  THE  RATIONAL  POLYNOMIAL 
APPROXIMATION  TO  THF  NCPMAL  PROBABILITY  FUNCTION. 
TA(X)=SOPT(ALOG(l.  0/X**2)) 

APX<X)=TA<X)-<2.515517f0.802953*TA  <XU0 . 01032 4*T  A  (X ) •  *21  / 
l<1.0*1.432?99*TA<X)+0.199269*TAm  **2*  0 . Q0130b*T  A  (XI **3) 

L0  =  MDSTR*-1  '  . 

GO  TO  (100,200,300)  ,IOI?TR 

100  IF (OME AN) 111, ill, 112 

111  DMEAN=0.407  '  ,  i 

S0=4, 0 

112  IFtNOSTP-il 101,101, 102 

101  PS<l)=OMEAN*l. 0E-6  , 

C5=S0**5  1 

DIAM<1)=DMEAN*C5 

0TAM<2)=0MEAN/C5  .  , 

FMASS (1) =1. 0 
GO  TO  400 

102  BARMU=ALOGIOMEAN)  1  l 

SIGMA=ALOG(SO) 

BARMU=BARMU*-T.  *SIGMA*»2  , 

FRAC=1.0/FLOAT<NDSTF)  , 

OO  103  NO=l,NOSTR 

103  FMASS (NO) =FRAC 

NH=NDSTR/2  < 

OO  104  1=1, NH  ‘ 

PRB=FLOAT(I)»FRAC 

OIAM(IH)=BARMU«-APX  <PRP)»SIGMA  , 

J=NOSTR-I*l  1  1 

104  OIAM<J)=BARMU-APX<PFB) »SIGMA  f 

FOR  THE  2  EXTREME  INTERVALS  THE  AVERAGE  DIAMETER  IS  ' 

ASSUMED  TO  »E  AT  HALF  A  MASS  FRACTION  FROM  ZERO  AND  ONE 

PR8=FRAC/2.0  ' 

PS(1)=BARMU*APX(PRR)*SIGMA 
PS(NOSTP)=BARMU-APX <PRBI*SIGMA 
OIAM(l)=2.*PS<l)-OIAM<2i 
0IAM<L0>=2.*PS<N0STF)-DIAM(MDSTR)! 

'  i 

CALCULATE  MEAN  DIAMETERS  FROM  BOUNDARY  VALUES. 

JsNOSTR-1 

IF<J-1)107,107,105 

105  DO  106  1=2, J  ■ 

106  PSm=0.5*(OIAM<mciAM<I+m 
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46 
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47 
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56 

OSTBN 

57 
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107 

00  108  1  =  1 , NDSTP 

OSTRN 

59 

OIAMfI»=FXP(OIAM<I) ) 

ORTON 

60 

100 

PS(I)=EXP<PS<I)1*1.0E-6 

OS  TON 

61 

DIAH(LO)=EXPCOIAM<LCn 

OS  TON 

62 

GO  TO  400 

OS  TON 

63 

200 

IF1EXP0-4. 01201, 202, 202 

OS  TON 

64 

202 

WPITE < ISOUT ,2001) 

OS  T«N 

6R 

i 

call  EXIT 

OSTRN 

66 

2001 

FORMATf*l*,lX,  *EXPONENT  OF  POWER  LAW  POWER  LAW  PARTICLE 

SIZE  FRE QUOSTBN 

67 

1ENCY  DISTRIBUTION  GT.  OP  EQ.  4.0*1 

OS  TON 

68 

201 

IF(NOSTR-11 203,204,204 

OS  TON 

6° 

203 

NOS.TR=10 

OSTRN 

70 

204 

AN=FLCAT (NOSTR1 

ostrn 

71 

FRAC=1.0/.AN 

OSTBN 

72 

OO  205  1=1 ,NOSTR 

OS  TON 

73 

20*5 

F1ASS(I)=FRAC 

OSTRN 

74 

POW=1.0/(4.0-EXPOI 

OSTBN 

75 

DMIN=l6.'0*SSAH*FRAC/CPOW*CAY»ONS*3.1415qE6)  >**POW 

OS  TON 

76 

00  206  IJ=l,NOSTR 

OSTRN 

77 

,  aj=floahij»-i.o 

OSTRN 

76 

’  206 

OTAM(IJ1=tAN-AJ1**PCH*nMTN 

■•OSTBN 

79 

PS<NPSTR)=ONIN*0.5**POW 

OSTRN 

RO 

DljAM(LO>=PS<NOSTR»  **2/0!AH(N0STRl 

OSTRN 

91 

NO=NOSTP-l 

ostrn 

82 

i 

00  207  I J=1 ,NO 

OSTBN 

33 

,207 

PS<IJ)=SQRT<OIAM(I  J)*OIAM(IJ«-l)  > 

OSTRN 

84 

00  208  IJ=l,LO 

DSTBN 

35 

208 

DIAM<IJ|=1.OE»6*0IAF(IJI 

OSTRN 

86 

GO  TO' 400 

OcTRN 

37 

300 

00  301  1=1 ,NDSTR 

OSTRN 

88 

301 

PS<I)  =  0.5MDIAM<I)+CIA*MH-111*1.0E-6 

OSTBN 

89 

400 

REfURN 

OSTRN 

90 

END 

OSTBN 

91 
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SUBROUTINE  SHWINP 

SHWNP 

1 

SHWNP 

2 

REAPS  IN  SHCT  TIME  WINE  OATA  ABOVE  GROUND  ZERO 

SHWNO 

3 

SHWND 

4 

COMMON  /SET1/ 

SHWNn 

5 

1C AY  ,  OETIO(12> 

, CIA ‘1(201) 

, DMEAN 

,DNS 

,EXPO 

, SHWNP 

6 

2FMASS(2001 ,  TOISTR 

, IFXEC 

, IRISE 

,  ISIN 

, ISOUT 

, SHWN° 

7 

3NOSTP  , PS ( 20  0) 

,SD 

,  S  S  AM 

,TME 

,  TMP1 

,  SHWNP 

8 

4TMP2  ,T2M 

, USOIL 

,VPR 

»w 

, HEIGHT 

, SHWNP 

Q 

5ZSCL  ,  NHOOO 

» ZV (200 ) 

, V  X (20  0 ) 

»VY (200) 

SHWNP 

If 

PEADdSIN,  DNHOOO 

SHWNP 

11 

IF(NHOPO>  100,  100,200 

SHWND 

12 

100  WRITE  < ISOUT ,51 

SHWNP 

13 

GO  TO  300 

SHWND 

14 

200  READ  (IS  IN,  2M  Z  V(  J)  , VX(J),VT(J»  ,J 

=1 ,NHODO) 

SHWND 

lr' 

W°ITE (ISOUT ,3INH000 

SHWNP 

16 

WR ITE ( ISOUT ,41  <  Z V  <  J 1 

,VX(J) ,VY(JI 

, J=l,NHODO> 

SHWNP 

17 

300  RETURN 

SHWND 

IP 

SHWNP 

1 R 

1  FORMA  T  ( 15) 

SHWNP 

20 

2  FORMAT (F12 « 3,  2F12.3) 

SHWNP 

21 

3  FORMAT ( Fi F , SX , FW INO 

HODOGRAPH  AT 

GROUND  ZERO*,tOX,*NHODO  =  F , 

I57/1SHWN0 

2? 

11 Y , FVECTPR  ALTITUDE, 

ZV( J) *,16X, 

FVX(J)F,24X 

,FVYC  J)F) 

SHWNP 

23 

4  F0RMAT(3(16X,E13.6)  ) 

SHWND 

24 

5  FORMAT (FIE, OX, FSHOT - 

TIMF  WTNOS  HAVE  NOT  BEEN  SPEOIFIFDF) 

SHWNP 

25 

END 

SHWNO 

26 
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